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LABORATORY TEST ON CYCLIC 


C.D.U. 550,41 :542.1/4:624.1431.25 


«CLOSED CIRCUIT» 


PERCOLATION OF ROCKS (”) 


RESUMO 


Simula-se um maciço rochoso e um túnelem rocha, sujet- 
tos a condições intermitentes de molhagem-secagem, 

Primeiro submete-se o maciço ou túnel a uma percolação 
provocada por aspersão de wm liquido a temperatura conhe- 
cida, Em seguida di-se a circulação de uma corrente de ar 
quente a temperatura estabelecida, Por Am há um periodo 
de repouso, Este ciclo com três fases repete-se, automútica- 
mente, 

Ensaiaram-se simulações de maciços rochosos de sienito 
nefelínico e de basulto oltvínico e um túnel em basalto oli- 
vênico. 

Comparam-se os resultados obtidos ao fim de ensaios de 
1000 ciclos. 


I—- LABORATORY TEST PROCEDURE 


This type of test was intented to simulate 
either a rock mass or a tunnel in a rock, exposed 
to intermittent moistering — drying conditions 
(seasonal regime, actually). 

Fig. 1 presents the diagram of the device 
designed for the test. The rock to be tested, a 
coarse-granulated rock, is subjected to percolation 
caused by sprinkling. 

The percolating liquid (at first distilled water), 
at a pre-established temperature which is main- 
tained during the whole test, is pumped into the 
sample-holder. The liquid penetrates the sample- 


Dv L. AIRES-BARROS & R. C. GRACA 
Laboratório de Mineralogia e Petrologia 
| S, 1. — Lisboa 


SYNOPSIS 


This laboratory test was intended to simulate either a 
rock mass or a tunnel in a rock, exposed to intermittent 
moistereng-dryinç conditions. 

lhe percolating liquid ata pre-establishedl temperature 
which is maintained during the whole test, às pumped into 
the sample-holder, A fter sprienklina, an authomatre cam 
system a'urts the circulatrion of a warm air curcentorer 
the sample. During the third stage — rest — the somple 18 
maintained under room pressure and temperature conditiona» 

As during the investigation, the various types of rocks 
were aubmitted to 1000 eyeles, we obtarmed many valuable 


data un the eftects of accelerate leachiuna, 


holder through a rose, so that the sample («rock 
mass») is subjected to sprinkling («rain»). The per- 
colating liquid is then carried back to the tub from 
which it had started, and will be used again in 
a later cycle. 


As we can see, it is a closed circuit type per- 
colation, as the water is only renewed when 
considered unsufficient to assure the perfect 
circulation in the pipes and through the sample. 


The diagram in Fig. 1 illustrates the circuit 
covered by the water and the position occupied 
by the various devices designed to accomplish 
percolation (water heating and temperature-contro] 


(*) This study was undertaken within a research project on «GEOCHEMISTRY OF ROCK EXPERIMENTAL 
ALTERATION», sponsored by the 1.A.C. of the National Education Ministry. 


Received May 16th, 1971. 
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Fig. 1 — Diagram of the closed circuit percolator 


system, water pump and volume controlling valve). 

After sprinkling, an authomatic cam system 
starts the circulation of a warm air current over 
the «rock mass». This warm air current reaches 
the sample at 50º €, During the third stage — rest- 


--A1s  TEMPORIZER 


The system is authomatically controlled by 
means of temporizers (see diagram on Fig. 2). 
The automation system consist of two tempo- 
rizing units. There is one temporizer which allows 
temporization from one second up to 48 hours 


e... ue a e a A rm A se ST 


TEMPORIZER 


Hig. 2 — Diagram of the electrical connections 


the «rock mass» is maintained under room 
pressure and temperature conditions. 

Each three-stage cycle is followed by another 
three-stage cycle-water percolation, warm air 
circulation and rest. The cycles are authomatically 
computed. 


q18 


maximum and one cam system temporizer. The 
first temporizer controls a pumping device which 
propels a small, corrosive-liquids-proof centrifuge 
pump for a pre-determined period of time. This 
pump was meant to start the circulation of the 
percolating liquid through. the sample-holder- 
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After percolation, this same temporizer sets in 
motion the cam system temporizer, which con- 
trols the air heating system, the warm air being 
produce by a ventilation system which contains 
an electrical resistance prepared to bring the air 
up to the desired temperature. This cam system 
temporizer controls also the rest stage, and oncf 
a cycle is finished (water circulation, rest) it 
sets the first temporizer in motion and a new 
cycle begins. 

As we have already said, the percolating 
liquid must be periodically renewed, because of 
josses caused by evaporation. The percolating 
liquid is renewed whenever it reaches a given 
level considered as unsufficient for the full 
leaching of the «rock mass». 

As during the investigation, the various types 
of rocks were submitted to 1000 cycles, we 
obtained many valuable data on the effects of 
accelerate leaching. The sample was analysed 
before the test (for definition of density, chemical 
analysis, micropetrographical characteristics), and 
again when the test was over. 

The next chapter will tell a bit more about the 
rocks which have been tested (mepheline syenite 
of Monchique and olivine basalt of the Lisbon 
region). 

This test, roughly described above, can be 
transformed into an artificial percolation through 
a «tunnel» in a rock. In this case, we used 44 mm 
diameter test tubes of rocks in which a «tunnel», 
with 13 mm diameter and 113 cm” side area, was 
drilled. Through this «tunnel» operated, in conti- 
nuous stages, percolation, warm air insufflation 
and rest. In a similar way, the percolating liquids 
were also systematically analysed. After the test 
(1000 cycles) we proceed to the mineralogic- 
-petrographical analysis of the walls Of the 
simulated «tunnel». 


H — RESULTS 
IH.1 Nepheline syenite 


Nepheline syenite is an undersaturated rock, 
sodic-alkali-silicates-rich (mepheline, albite-oli- 
goclase and aegirine). A marked mobilization of 
the alkalies, especially soda was therefore to be 
expected. 

From Table 1, Fig. 3 was worked out, and the 
following conclusions derived : 
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a) The proportion of mobilized Na,O was 
higher than the proportions of any other leached 
material. This Na,O was produced either by 
nepheline, or by albite, or else by aegirine. 

b) We must stress that, in heeping with the 
mineralogical and chemical composition of the 
rock, the proportion of mobilized Na,O was 
about three times higher than the proportion of 
mobilized Ca O, and four times higher than the 
proportion of mobilized K,O. (Table 1). 


TABLE 1 
Nepheline syenite 
Granulate 
SiOs | Cao Naj0 | Ko0 Fe,0, 
Epa º 0><10"|9,><10"|9 >< 10" |" >< 10º| 1,9 >< 10º 
25 — 0,279 3420 | 1,125 | — 
29 — 0,315 3,627 | 1,877 | — 
52 — 0,688 3,906 | 2,629 — 
77 — o,g6o 4,321 2,700 — 
ITI | 0,322 2,149 5,580 | 2,500 0,064 
161 0,044 2,708 6,011 | 3,021 0,099 
178 | 0,931 3.115 6,362 | 3,031 0,145 
373 | 1,253 | 3,258 | 8097 | agro | o,t45 
soo | 1,253 | 4,835 13,695 | 3,783 | 9,790 
655 | 1,360 | 5435 | 14,651 | 3914 | 0,790 
8oo | r,360 | 5.518 | 15,792 4155 0,790 
845 | 1,360 | 5,963 17,129 | 4308 | 1,159 
1000 1,360 | 6374 | 17,767 | 4,328 | 1,521 


c) The high proportion of mobilized Ca O was 
due, not only to feldspars, but also to sphene 
or dolomite (of secondary genesis), also present 
in the rock. 

d) K,O was produced by sodic-potassium 
feldspars. 

e) The mobilized iron was produced by the 
mafic minerals present in the rock (aegirine, tita- 
nite and magnetite). 

t) While the intersticial cations migrated easily 
(Na' = Ca'+ a Fe*” ; Mg” * ), silica remained 
in the structural lattice of the silicates, and the- 
refore its mobilization rate was low. 


WH.2 Olivine basalt 
Olivine basalt is also an undersaturated rock, 


but a calk-alkaline one, rich in mafic minerals 
(augite and olivine) and also magnetite-rich. For 
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NEPHELINS SYENITE OLIVINE BASALT 


CRANULATE 


Fig. 3 


olivine basalt the investigation was made on gra- 
nulated rocks, and into a cylindric «tunnel». The- 
refore, it will be most useful to examine Table II 
together with Table III. 


TABLE 1 
Olivine basal! Fig. 4 
Granulate TABLE HI 
Olivine basalt 
| Sio, | Cao Na,0 | K 50 Fe,0, Cylindric tunnel 
Cycles 1103 95 10º %><10" "510" ar 0! nação! 10º Proportions as compared to percolation surface 

| SiO, | Cao Na,0 | K50 | Fes0 
122 | 0,2 0,6 51 | I,2 LI Cycles |, ' Jo J ; 3). e e 
09 | 36 [aulas lar | se) | |uxi0!tne1o!] My ><10?] > 10º] voe ao? 
285 45 | 1,6 8,9 2,5 59 21 1,8 2,2 08 | 04 0,9 
363 5.0 39 | Tor 32 7,7 87 r,8 3.3 19 1,2 45 
463 | 5.1 5,6 12,5 4,2 10,6 140 1,8 2,6 2,4 1,5 5,6 
516 5.7 | 6,1 13,2 4,9 12,2 174 3,8 II,O 86 | 36 10,9 
591 6,4 6,4 14,4 594 14,7 310 58 | 13,9 1,5 4,8 12,5 
652 6,5 75 15,0 5,8 15,7 438 7,3 25,1 15.8 6,7 15,2 
717 6,9 9,3 16,6 6,2 20,9 566 7,0 32,2 19,8 8,6 20,0 
775 75 10.5 16,9 6,4 21,5 665 9,2 35,8 20,9 9,0 23,2 
833 7,6 IL,5 17,4 6,5 23,5 729 [5,0 37,8 22,0 | 9,7 24,4 
gr 7,8 13,3 17,8 6,7 23,5 789 16,4 40,5 22,7 IO,1 25,6 
041 9,2 16,0 18,2 6,9 25,3 841 16,8 45.0 23,1 10,5 27,7 
Iooo | 94 17,9 18,5 71 32,3 1000 17,1 499 24,3 11,7 29,8 
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CYLINDER 


| ATTTA* 
e — [HÃO 
Pena pin f | / | 

SO p—— des: Os Tr — e 7 | | | / ] = = 


a+* 
..* 


1000 
CYCLES 


Fig. 5 


a) Notice that leaching was higher for the 
percolation accomplished in the «tunnel» than 
for the percolation accomplished on the granulate. 
Percolation in the «tunnel» was accomplished 
under pressure, causing the removal of larger 
quantities of material, and this may explain the 
fact. Besides, despite the 16 mm thickness of the 
«tunnel», the temperature around its external 
wall, when the percolation (both liquid and air) 
took place, reached 10 to 15º € over room tem- 
perature. Moreover, the external walls of the 
«tunnel» presented occasional fine liquid exuda- 
tion. 

b) Apart from the figure for Fe, O,, which was 
similar in both «tunnel» and granulate, the 
leaching figures for the «tunnel» were between 
about 1.5 and 3 higher than for the granulate. 

c) In general, Ca O leaching is higher than 
K,O and Na, O leaching; we think this happens 
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because the rock is calk-sodic-plagioclasses-rich. 
For the same reason the figure for mobilized 
Na, O was twice the figure for leached K, O. 

d) The quantity of Fe, O, supplied is unders- 
tandable, as olivine basalt is rich in mafic 
minerals. 

The findings for olivine basalt prove that the 
percolation was much stronger in the cylindric 
«tunnel» than on the granulate. To this has 
surely contributed the persistant artificial leaching 
in the «tunnel», in which the walls «sweat», 
denouncing a deep penetration of the percolating 
liquid. 

In general, the figures concerning the mobi- 
lized material in the «tunnel» were about two or 
three times higher than on the granulate. The 
only exception to this rule was Fe, O,. It happens 
that basalt is very rich in magnetite microlites 
which are homogeneously distributed all over 
the rock mesostasis. That explains why the 
figures for leached Fe, O, were similar on both 
occasions. On the contrary, feldspars, suppliers 
of CaO, Na,O and K,0, are irregulary distri- 
buted. These minerals may occur either in micro- 
lites or in phenocrysts, and the phenocrysts 
expose to leaching a smaller surface per volume 
unit, and this explains its low removal as the 
water penetrates as softly as rain in the granu- 
late; the same cannot be said for the artificial 
penetration in the «tunnel». 

In short, whenever a mineral occurs in a dif- 
fused and regular distribution and in small crys- 
talls, exposing to leaching a large surface per 
volume unit (usually in fine-grained rocks), it 
looses larger quantities of its contents than the 
same mineral occuring in large cristals and irregu- 
lar distribution. In the granulate, the liquid goes 
softly through, thus leaching the microlites in 
the perphery of the fragments, but it does not 
reach the core of the fragments. In the «tunnel», 
with «sweating» external walls, the liquid, at 
a high temperature, goes deep into the fragment, 
being therefore of an extreme agressivity. 


I.3 Comparison between nepheline syenite 
and olivine basalt 


To compare the behaviour of these two gra- 
nulated rocks, we simulated rock mass under 
warm «rain» conditions alterning with intense 
dryness conditions, and Table IV was worked 
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out. This Table indicates not only the main oxi- 
des in the sound rocks, but also the proportion 
of leached oxids after 1000 cycles. 


TABLE IV 


Olivine basalt 


Nepheline syenite 
(Lisbon) 


(Monchique) 


Leached pro- 
jduct (1000 cyc.) 
5 


ti | Leached pro- 
'duct (1000 cyc.) 


Sound rock 


o rock "y 


; | 09 >< 10 | 9/y >€ 10º 
Si O, | 54,36 | 2,5 | 4213 22,3 
Fe, Os | 2.57 59,2 12,13 266,3 
Ca O NH 306,3 10,57 | 169,3 
Na, O | 9,02 197,0 | 2,29 | 807,9 
K, O 7,20 60,1 | 0,77 | 922,1 


From Table IV some conclusions may be 
derived: 


a) Olivine basalt, which has a lower content 
of silica as compared to nepheline syenite, exports 
nine times the silica exported by nepheline 
syenite, 


b) Olivine basalt, poorer in alcalies than nephe- 
line syenite, exports about four times more Na, O 
and fifteen times more K, O than the latter. 
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c) Though olivine basalt has higher contents of 
CaO than nepheline syenite, this one exports 
about twice the lime exported by olivine basalt. 


d) The total sum for iron mobilized by olivine 
basalt is five times higher than the sum exported 
by nepheline syenite. 

The explanation for these facts apparently 
discordant, is that the feldsparts in the basalt 
are more basic and aluminous than the felds- 
parts in nepheline syenite, and therefore liberate 
more easily and readily their intersticial cations 
(Na", K', Cat Thus, the leaching of the 
cations disturbs the equilibrium of the structural 
network, which collapses simultaneously libera- 
ting silicon. 

It is therefore proved that the basic rocks, 
basalt for instance, mobilize a larger quantity of 
intersticial cations or even sillica as compared to 
alkaline rocks, which are richer both in alkalies 
and silica. 

The figures for CaO are different, as nepheline 
syenite is a better supplier of lime than basalt. 
The petrographical analysis shows that nepheline 
syenite contains secondary dolomite. This explains 
the high figure for CaO exported by nepheline 
syenite. 

Being rich in mafic minerals, olivine basalt 
was, Of course, a better supplier of iron as com- 
pared to nepheline syenite, which has poor 
contents of these minerals. 
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NOTA SOBRE A APLICAÇÃO DAS TRANSFORMADAS 
DE FOURIER E DE LAPLACE AOS CIRCUITOS” 


RESUMO 


Pretende este artigo esclarecer algumas dúvidas ausci- 
tadas pelos estudantes quanto à aplicação aos circuitos das 
transformadas de Laplace e Fourier, 


1— INTRODUÇÃO 


Entre os estudantes a quem foram apresenta- 
das as Transformadas de Fourier e de Laplace é 
frequente aparecer a seguinte dúvida: na aplica- 
ção destas transformadas a circuitos, em que 
circunstâncias é mais aconselhável a utilização 
duma ou de outra; a semelhança entre as duas 
transformadas justifica esta hesitação. Na pre- 
sente nota, a publicar numa revista de estudantes, 
faz-se uma revisão do assunto, tentando respon- 
der, pelo menos parcialmente, àquela interroga- 
ção. 


2— A TRANSFORMADA DE FOURIER E A 
SUA INTERPRETAÇÃO EM TERMOS 
DO ESPAÇO TEMPO-FREQUÊNCIA 


Vai-se usar aqui a seguinte definição de Trans- 


formada de Fourier (TF) 


0) TELEQ]=Fi)= DD [e ita 


— OO 


Do ponto de vista dos circuitos f(t) represen- 
tará a descrição ao longo do tempo duma gran- 
deza intensiva mensurável e será portanto uma 
função real da variável real t. Nestas condições 
e se a idealização do problema admitir que 


oo 


(2) | 1£(9| dt <<os 
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SYNOPSIS 


This article is intented to clear up some doubts obout the 
cireuitapplications of the Laplace and Fourier transforma, 


então F(») é uma função (complexa) limitada e 
contínua sobre todo o eixo real da variável w e 
[F(»)| - O quando |»| — co/1]. Juntando ainda o 
facto de que por ser f(t) real se tem F(— w) = 
==F*(m») vê-se que a TF permite passar duma 
representação no tempo, com t G (— co, + «), para 
uma representação na frequência angular » = 27f 
em que a frequência f é uma grandeza real e 
positiva. 

Por outro lado, se a função f(t) acima consi- 
derada for ainda uma função seccionalmente 
contínua, com um número finito de discontinui- 
dades de 1.º espécie em cada intervalo finito, 
mostra-se [1] que se lhe pode dar a seguinte 
representação em cada ponto t: 


(3) = Etr +tt-o]= [| Fe)el'tdo 


em que F(») é a TF definida em (1) e 


f(a+o) =lim f(t) 


toa[!Z: 
t<a 


A expressão (3) pode ser considerada a ope- 
ração inversa de (1) e pode além disso servir de 
base para a seguinte interpretação: em cada 
instante t (excepto eventualmente nas suas 
discontinuidades) a função f(t) pode ser sinte- 
tizada à custa duma sobreposição de sinusoides 
de frequência w e amplitude F(») dw em que 
F(») é uma densidade (complexa); estende-se 
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assim a funções não periódicas um resultado já 
conhecido para as funções periódicas (Série de 
Fourier). Determinar F(») por (1) corresponderá 
portanto a fazer a análise espectral de f(t) 
enquanto que por (3) se faz a síntese da função 
a partir da sua densidade espectral. Note-se que 
no presente contexto não têm representação 
espectral certas funções que correspondem a 
idealizações simples e naturais; é o caso da fun- 
ção de Heaviside (representando o estabelecimento 
duma corrente contínua) e o das funções da 
forma f(t + T)=-f(t) para qualquer valor de 
tC (0,0) —- em ambos os casos não se verifica 
a condição (2), sendo |f(t)|F0 (mas limitado) 
quando t —» co (*). 


3— APLICAÇÃO DA TF AOS CIRCUITOS 


Passa-se agora à aplicação da TF aos circuitos: 
só se consideram aqui circuitos lineares (para os 
quais é válido o princípio da sobreposição) e 
invariantes no tempo (invariantes sob uma trans- 
lação em t). Nestas condições pode-se simplificar 
muito esta discussão sem diminuir a generalidade 
das conclusões considerando apenas quadripolos 
desprovidos de qualquer fonte independente, em 
que um par de terminais (saída) está numa situação 
de carga passiva (isto é, que até qualquer instante t 
absorveu pelo menos tanta energia quanto a que 
devolveu) mas arbitrária, e em que se liga ao 
outro par de terminais (entrada) uma fonte 
independente; a aplicação sucessiva do teorema 
de Thévenin (ou do de Norton) permite reduzir 
a esta situação qualquer circuito linear com 
várias fontes independentes, enquanto que a 
sobreposição permite sintetizar em qualquer dos 
ramos do circuito o efeito das várias fontes 
conhecido o efeito nesse ramo de cada uma delas. 
Ao conjunto quadripolo + carga dá-se aqui a 
designação de sistema. 

Tendo estabelecido esta situação, pode-se agora 
fazer uma experiência: tome-se para fonte um 


(*) Por vezes atribui-se uma descrição espectral a estas 
funções, calculada por meio duma passagem ao limite 
a partir de funções obedecendo à condição (2); os espec- 
tros assim obtidos exibem polos sobre o eixo real e a sua 
inversão deve ser feita intentando o caminho de integração 
usado em (3) de modo a deixar os polos no semi-plano 
superior no plano « complexo — o que equivale a usar um 
caminho de integração como o descrito no parágrafo 4, 
com c >0. 
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gerador sinusoidal, isto é, impondo na entrada 
uma grandeza intensiva (tensão ou corrente) dada 
por 


aç=Açelt! ,g(-co 465) 


em que 


A: constante complexa (amplitude complexa) 
w==2 f: número real positivo (frequência angular) 


Sendo o circuito linear e invariante no tempo 
observa-se que qualquer das grandezas intensivas 
na carga exibe aquela mesma forma de variação 
no tempo, com a mesma frequência mas de 
amplitude complexa em geral diferente ou seja 


RA jwt 
aç=A e 


Entre as duas amplitudes complexas, A, da 
excitação e A, da resposta existe uma relação 
imposta pelo sistema: 


(4) T(u) = ÀS 
A 


e 


tendo-se explicitado que essa relação é em geral 
função da frequência. Por outro lado está impli- 
cito nesta definição que t€ (— co, + o9) atra- 
vés da definição da fonte, o que implica que 
o circuito manteve permanentemente a sua estru- 
tura. 

Pode-se agora utilizar a interpretação dada 
a (3) em termos da síntese da função por meio 
de funções sinusoidais para obter a resposta 
f.(t) do sistema a uma excitação imposta f, (t). 
Assim ter-se-á 


86) 


(5) fo(t) = ) T(w) E. (m) elet do 


o 
em que 


E. (0) = TE [E (9] 


O integral (5) existe se T(») for limitado, ou 
seja, |T(m)|<M, visto que sendo assim se tem 


| DO 
co M f A (19) e! nei d (1)| === 


— (O) 


fre F ()e!“' das 


— (CX) 


= M | fo(t)| 
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Admita-se agora que a resposta g (t) do sistema 
à função de Dirac satisfaz à condição 


(6) fiswia<es 


Existe então uma decomposição espectral de 


eg (t) em sinusoides de frequências reais. Para 
a obter recorre-se a (5) sabendo que a TF [9 (t)] = 
==1'27; vem 


E ()= [ riu) "E da 
2 


e por comparação com (3) obtém-se 
1 

TELg (t)] = T (0) 
2% 


A condição (6) garante que T (») é limitado 
e portanto a existência de (5) para qualquer 
excitação f. (t). 

Tanto g(t) como T(m) — que se vai designar 
por função de transferência (FT) — são carac- 
terísticas do sistema (quadripolo + carga). 

Pode-se agora concluir sobre a aplicação da 
TFa circuitos: parte-se dum sistema linear e inva- 
riante no tempo, imutável em tC(— co, + 00) 
e ainda para o qual se verifica a condição (6). 
Está-se então em condições de definir uma FT 
que se determina pelo cociente das amplitudes 
complexas de sinais sinusoidais; essa FT é con- 
tinua e limitada sobre o eixo w real. Formou-se 
deste modo do ponto de vista teórico uma base 
para o estudo das transformações sofridas pelos 
sinais através de sistemas em que se utilizam as 
características de frequência duns e doutros 
— um exemplo típico é o dos filtros cuja FT é 
dimensionada de modo a aceitar ou a rejeitar 
determinados sinais definidos pelo seu espectro. 
Por outro lado a própria definição (4) serve de 
base para a determinação experimental da FT: 
o sistema é excitado por meio dum gerador dando 
origem a um sinal «sinusoidal» cuja frequência 
é definida tanto mais precisamente quanto mais 
longo for o tempo de observação — trata-se na 
realidade dum procedimento universalmente 
adoptado. 
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4 — TRANSFORMADA DE FOURIER UNILA- 
TERAL 


Passar-se-ia agora a considerar a transformada 
de Laplace (TL), definida por 


co 


TL[E(t)] es | f(t)eT Pt dt 


É no entanto mais conveniente aqui lidar com 
a transformada de Fourier unilateral (TEU) defi- 
nida por 


(7) TEU [Ely] == E | E) eiot ar 


e que difere portanto da TL pela constante mul- 
tiplicativa 1/27 e por se ter feito p=)“; com 
esta mudança não se altera as conclusões a que 
se pretende chegar e por outro lado facilita-se 
a comparação com a IF. 

Acentua-se em 1.º lugar que E” (») é deter- 
minado pelo conhecimento de f(t) em tG (0,00); 
este facto implica o conhecimento de valores ini- 
ciais para o cálculo da TFU da derivada e do 
integral de f (t). 

Por outro lado a condição (2) de integrabili- 
dade absoluta de f(t) é relaxada, passando-se 
só a pedir a integrabilidade absoluta de f (t) em 
qualquer intervalo finito e que seja 


t(t)=0 (e!) 


(a > 0) quanto t— 9; F ” (w) é então analí- 
tica no semi-plano [2] 


Img [»] <—a 


e é limitada e contínua sobre a linha Img [1] = — a. 

Para inverter esta transformada recorre-se a 
uma função u (t) = f(t) et, com c>a, ten- 
do-se arbitrado que f(t)=-0 em t < O; é agora 
aplicável a esta nova função o integral de inver- 
são (5), vindo então 


O — 1€ 
. 


D  dy=] 


— E — |O 


ET (1) elit do 


em que a condição c >> a garante que o caminho 
de integração existe na região de analiticidade 


de ET (1). 
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É evidente que nestas condições se perde 
a interpretação dada acima à TF de consistir na 
análise de f(t) em funções sinusoidais de fre- 
quência real. Em compensação alargou-se a classe 
de funções para as quais é aplicável este tipo 
de transformada; é de esperar que também se 
tenha alargado a classe de sistemas (no sentido 
aqui usado) que a transformada pode tratar. 


5 — APLICAÇÃO DA TFU AOS CIRCUITOS 


Na aplicação da TFU a um sistema (qua- 
dripolo + carga) linear admite-se que cada uma 
das grandezas intensivas intervenientes no sis- 
tema tem TFU; utilizando então a linearidade da 
transformada obtém-se, por eliminação das outras 
grandezas que não interessam ao problema, 
uma expressão a que se pode dar a forma 


(9) EO (0) =0(0)FÇ (1) + 0 (1) (to) 


em que F/ (») é a TFU da função de excitação 
Fo (t), Fl (») é a TFU admitida para a resposta 
F.(t), O (») é uma função do sistema e I (1) é uma 
função linear das condições iniciais que surgem 
ao calcular-se as TFU das derivadas e integrais 
das várias grandezas intensivas no quadripolo ; 
pode-se mostrar que esta última função é nula 
no caso do sistema estar em «repouso», isto é, 
na ausência de actividade energética no sistema 
anterior à aplicação do sinal de excitação (impli- 
cando valores iniciais nulos) — note-se que esta 
situação em alguns casos pode ser irrealizável na 
prática em virtude da inevitável presença de 
fontes de ruido que tenderão a manter a activi- 
dade energética independentemente de qualquer 
excitação exterior. 


Fazendo agora uso de (8) vem 


«x — 1 € 
(10) E, (= | [90) FÉ (1) +90) Eo)Je!"t dio 
— vo — jé 


em que « é escolhido de modo que o percurso de 
integração passe abaixo de todas as singulari- 
dades da função integranda. 

Para poder interpretar “ (w) vai-se considerar o 
caso em que o sistema está em repouso e por- 
tanto I(»)==0 e admite-se que a excitação é 
feita por meio duma função de Dirac no instante 
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; 1 E 
t,; neste caso F; (9) = “elo, E (t) = g (t), 
2% 
vindo 
co ic 
E 1 tu ' Lts 
Blt—to)= | 90) elo et lot dy 
J 2% 
-— o JC 


Admitindo ainda que a função do sistema 
0 (») não tem singularidades em Img [»]<.O e 
que portanto se pode fazer c==0, a comparação 
com idêntica situação para a TF mostra que 
nestas condições 9 (m) == T(»). 

No entanto do ponto de vista da TFU é admis- 
sível que !(») tenha singularidades em Img [1] 
<. O (a distância finita) e que g (t) possa não ser 
absolutamente integrável em todo o intervalo 
(0,); por exemplo, a existência dum par de 
polos situado em v= +», —ja (a >0) conduz 
a uma resposta g (t) da forma 


k e” cosm,t = constante 


que óbviamente não satisfaz a (6) —-mas que é 
de ordem exponencial quando t —» «=. É de resto 
possível demonstrar que a resposta dum sistema 
a uma excitação de amplitude limitada ( | fu (t) 
< M com M finito) é também de amplitude 
limitada se e só se for observada a condição (6) 
para g (t) ( 3], pag. 670); nestas condições diz-se 
que o sistema é «estável» pelo critério acabado 
de enunciar. 

Verifica-se assim que a TF tal como foi consi- 
derada nos parágrafos anteriores terá a sua apli- 
cação restringida a sistemas lineares estáveis 
enquanto que a TFU é ainda aplicável a sistemas 
que poderão não ser estáveis. 

Finalmente, como consequência imediata do 
facto de só se considerar t > O na definição da 
TFU, observa-se que esta transformada está par- 
ticularmente bem adaptada ao caso em que se 
pretende determinar a resposta dum circuito que 
de algum modo é bruscamente alterado no ins- 
tante t=-0 — dão-se dois exemplos na figura 
junta. Para a aplicação da TFU só interessa a 
forma que o circuito assume para t >0, expri- 
mindo-se como condições iniciais o estado em 
t==+ O resultante da energia que, por um meca- 
nismo que não é necessário especificar, se tenha 
acumulado previamente (em t 0) nos elementos 
reactivos. Uma situação destas pode ser estu- 
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dada utilizando a TF desde que se suponha que 
o circuito se manteve desde t == — co com a forma 
que assumiu em t== + 0 e se especifique as fon- 


a SD ad 
Lo TLT TT. 


LET 


t=-0 t=+0 


Fig. 1 


tes de modo a ter f.(t)==0 para t<0-—nos 
exemplos da figura a bateria seria substituída 
por uma fonte dando origem a uma função de 
Heaviside (*). No entanto seria necessário que O 
circuito estivesse inicialmente em repouso ou 
então que fosse especificado o mecanismo pelo 
qual a situação observada em t==+ O tinha 
sido atingida. 


6 — REALIZABILIDADE FÍSICA 


Faz-se agora referência a um critério mencio- 
nado muitas vezes quando é tratada a aplicação 
da TF a circuitos. 

Considera-se então uma função f(t) que se 
admite satisfazer a (2) e prolonga-se a sua TF 
para fora do eixo w» real; para isso faz-se: 


E (0) =" (m) + E” (6) 


com 


oo 


E" ldem [ Eq gr 
= 


O 


E” (ij + | E(t) eta 
2: 


— (O) 


Desde que f(t) satisfaça a (2) será [2] 


(*) Ver a nota da página 424, 
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E" (») e F— (») contínuo e limitado sobre o eixo 
real 

E” (w) analítico em Img[»] SO e |F*(w)|- 0 
quando ||» cc em Img [»]<0 

F (») analítico em Img [»]>0 e |F  (v)/-0 
quando |» |» co em Img [»] >0 


Mas se 


É (t) =0 


em t0 


vem 
EF (19) = F ” (10) 


e portanto será F(») analítico em Img [»] O 

além de contínuo e limitado sobre o eixo real. 
Aplica-se este resultado a um sistema que 

satisfaça a (6), relembrando que é então 


+ TW)= TE] 


2 x 


Verifica-se que se para este sistema for g (t) = 
= 0 emt< 0 será T(») analítico em Img [»] < 0. 
Reciprocamente, se T(wm) for analítico em 
Img [»] JO e ainda |/T(»)| » O quando 
será g (t)= O para t< 0, como o mostra a apli- 
cação do teorema de Cauchy ao contorno repre- 
sentado na figura junta. 


| —+ OO 


Fig. 2 


Este resultado caracteriza o sistema, visto que 
envolve as suas funções características T(m») e 
g (t) e é tido usualmente por traduzir a «realiza- 
bilidade física» do sistema no sentido de que o 
sistema satisfaz ao princípio da causalidade: a 
resposta g (t) à função de Dirac aplicada em 
t==0 só toma valor diferente de zero em t > 0. 
O mesmo é ainda verdade para a resposta a 
qualquer outra forma de exitação iniciando-se em 
t==0 — como se infere de (5). 
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A ausência de singularidades de T(m») em 
Img [»] < O constitui portanto um teste elementar 
de realizabilidade física que eventualmente se 
pode tornar útil em problemas de síntese de 
circuitos a partir de propriedades desejáveis : um 
exemplo quasi trivial é o da inversão de filtros 
que não é realizável se, como é evidente por este 
critério, a FT do filtro a inverter exibir zeros em 
Img [»] <<. 

Acentua-se no entanto que este resultado foi 
obtido tendo em mente sistemas satisfazendo 
a (6), pelo que se excluem deste critério os 
sistemas instáveis (que precisamente não verificam 
(6), como se mencionou no $ 5), Na realidade 
estes últimos sistemas são fisicamente realizáveis 
no sentido de satisfazerem à causalidade, visto 
que só iniciam a sua actividade após a acção 
duma percussão inicial, ainda que depois a sua 
evolução acabe por se tornar independente dessa 
percussão. As próprias condições iniciais podem 
fazer o papel da percussão inicial, como é evidente 
de (10), mas este caso ainda pode ser interpretado 
como não violando a causalidade. Por outro lado, 
não é possivel estabelecer por meio da TFU um 
critério de realizabilidade física como se fez com 
base na TF visto que a escolha da constante « 
em (10) é feita de modo a deixar em Img [»| = —c 
todas as singularidades do integrando. 

Convém ainda notar que se pode mostrar para 
os sistemas de parâmetros concentrados em que 
as singularidades de º (w) se limitam a polos que 
a instabilidade está ligada à noção de quadripolo 
activo: num quadripolo activo se forem (v,,i.) 
as grandezas intensivas de entrada (positivas 
para o fluxo de energia entrando no quadripolo) 
e (v.,i,) as grandezas intensivas de saída 
(positivas para o fluxo de energia saindo do 
quadripolo), observa-se no caso presente que 
existe um instante t, tal que para t > t, se tem 


AD | IveOiO-v (Di ()dr<o 

(ou seja, a partir de t, o quadripolo cede à carga 
uma energia superior à que recebeu da excitação 
— evidentemente à custa da transformação da 
energia cedida por uma fonte exterior); para 
verificar (11) basta admitir que a instabilidade 


foi iniciada por uma perturbação de duração 
t 
finita pelo que El vo i. dY assume um valor cons- 


— OO 
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tante a partir dum dado instante, enquanto que 
t 


se O (w) tem polos em Img [»] “O fvs i; dz 

— (CH) 
cresce indefinidamente visto que esses polos 
conduzem a soluções em que intervém o factor 
et com a >0 (*). 

Conclui-se que o critério de realizabilidade 
física sugerido pela TF permite reconhecer um 
sistema fisicamente realizável (no sentido de satis- 
fazer à causalidade) mas não permite distinguir 
entre um sistema que não é fisicamente realizá- 
vel e um sistema fisicamente realizável mas 
instável (**), a não ser que se saiba de 
antemão se a resposta do sistema à função de Dirac 
satisfaz ou não a (6)—-mas este conhecimento 
pressupõe normalmente que g(t) é conhecido, 
o que permite classificar imediatamente o sistema 
do ponto de vista da causalidade, tornando inú- 
til a aplicação do critério. 


7 — CONCLUSÃO 


Da análise feita nos parágrafos anteriores veri- 
fica-se que, do ponto de vista da sua aplicação 
ao estudo dos circuitos, é difícil estabelecer uma 
distinção acentuada entre a TF e a TEU (ou, o que 
é o mesmo, a TL). 

Observou-se no entanto que a TFU é aplicável 
a uma classe mais geral tanto de circuitos — visto 
que inclui os circuitos instáveis em que uma per- 
turbação limitada inicia uma resposta cujo cres- 
cimento é ilimitado — como de excitações — que 
podem ser representadas por funções de ordem 
exponencial quando t —» co, Esta maior genera- 
lidade justifica a predominância adquirida pela 
TL na bibliografia moderna sobre circuitos — como 
exemplos menciona-se [3] e [4]. 

No entanto a TF presta-se a uma interpreta- 
ção em termos de análise e de síntese de sinais 
com base em sinusoides de frequência real; esta 


(*) Há que admitir que a carga é um dípolo passivo 
onde há dissipação de energia pelo que não só é 
t 


[v. i; dz > 0 como também se pode mostrar que v, ei, 


+ 
— > 


acabam por assumir o mesmo tipo de variação [4]. 


(**, Tem interesse mencionar neste ponto que na defi- 
nição de realizabilidade física dada em [5' (pág. 108) se 
exige que não existam modos de oscilação cujas ampli- 
tudes cresçam indefinidamente com o tempo. 
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interpretação é largamente utilizada do ponto de 
vista teórico — na análise e projecto de sistemas 
no que diz respeito à ocupação de bandas de fre- 
quência (o que é particularmente importante em 
comunicações) e correspondentemente na deter- 
minação experimental das características desses 
sistemas. Porém a coincidência da função do sis- 
tema definida pela TFU (função 9 ('»)) com a fun- 
ção de transferência (T (»)) precisamente para os 
sistemas a que a TF é aplicável, implica nova- 
mente que, mesmo deste ponto de vista, é indi- 
ferente pensar em termos duma ou de outra das 
transformadas. 

Num único ponto a TF se distingue da TL: 
trata-se do critério de realizabilidade física — que 
não é acessível à TL visto que esta transformada 
só é definida para t =0 e ao construir-se o seu 
integral de inversão (8) se impõe f(t) = O para 
t<* 0. No entanto, como se viu, este critério é de 
aplicação limitada ainda que se possa tornar útil, 
em particular quando se queira rejeitar quadri- 
polos irrealizáveis ou instáveis. 

Para terminar observa-se que se cingiu estas 
conclusões ao quadro das definições dadas nesta 


nota, mas que correspondem de resto ao campo 
de aplicação usual daquelas transformadas no 
domínio dos sinais determinísticos. 
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V— INTRODUCTION TO 
LASER DOPPLER — SHIFT ANEMOMETRY (*) 


RESUMO 


Resumem-se os mélodos ópticos propostos para o estudo 
do escoumento de fluidos. Indicam-se os seus fundamentos, 
vantagens e desvantagens. 

Em particular desenvolve-se a teoria do anemômetro de 
alusser», baseado no efeito de Doppler. E, descrevem-se os 
metodos de processamento des dados, em estudo no Imperial 
College. 


INTRODUCTION 


Turbulence is a phenomenon which is well 
known in engineering practice, since it is asso. 
ciated with many problems engineers have to 
overcome. The awareness of the desirable and 
undesirable effects of turbulence is still growing. 
It is generally recognised, however, that future 
progress in understanding and controlling many 
technical processes associated with turbulent 
flows may result from information regarding the 
local mean velocity and the fluctuating compo- 
nents. This has increased the incentive to develop 
suitable measuring techniques to obtain the 
required information. 

To date, experimental investigations of turbu- 
lent flows have been restricted, almost exclusively, 
to low-speed subsonic boundary layers. Flows of 
this kind have been extensively investigated by 
means of hot-wire anemometer techniques and 


(*) Received May 20th, 1971, 
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SYNOPSIS 


A short summary of optical measuring techniques, 
proposed for fluid flow investigations, is given, The basic 
principles of operations are explained for the different 
methods and some advantages and disadvantages are 
pointed out. 

The main part of this paper deals with laser Doppler- 
shift anemometry aqvhich is a new and very promising 
meth d for local optical velocity measurements, With this 
technique the Jrequency shift of light, caused by scattering 
from small partieles moving with the fluid, is measured, 
This acattered light shows an appuarent change in frequency, 
the Doppler-shift, which can be utilised for velocity measu- 
rementa, 

The problems that oceur when an optimal optical system 
is designed are cutlined. Four different methods of data 
processing under investigation at Imperial College are 
described, For gas flow investigations, the range of suitable 
tracer particles is given and a table is presented to allow 
the solution of suitable ascaitering targeis for other fluids, 


a considerable volume of experimental data has 
become available. However, while hot-wire ane- 
mometry has proved to be an invaluable method 
for such flow cases, in fact virtually the only 
one to find general adoption, its application has 
not been generally extended to separated flows, 
high-speed flows, turbulent diffusion flames, etc. 
Experimental investigations of such flow situations 
are of current interest. However, they are inac- 
cessible to hot wire anemometry due to the pre- 
sence of shock waves, probe interference, strong 
fluid property gradients, high turbulence fluctua- 
tions, etc., resulting in difficulties of signal inter- 
pretation. 

The above problems, and the fact that probe- 
-type instruments cannot be used in extreme 
temperature environments, suggest the use of 
optical velocity measuring techniques. Such 
techniques are well known and it seems appro- 
priate to give a short summary of these methods 
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and their applications before describing the basic 
principles of Doppler-shift anemometry. The 
latter technique is a new and very promising 
approach to obtain local velocity measurements 
in laminar and turbulent flows of transparent 
fluids. This new method has great advantages in 
comparison with other methods, see (1). The 
most important ones are given within this paper, 


. Pa E 
| | | Beam 
=) | ' [TER 


fc — INTERFEROMETEY 


Fig. 1 — Sketches of optical arrangements of standard 
optical velocity measuring techniques 


Fig. 1 shows the basic optical arrangements 
of three methods that have been employed in 
the past for optical velocity measurements. These 
methods, Shadowgraphs, Schlieren methods, and 
Interferometry, are standard techniques in fluid 
flow research. All of them are well described in 
in many textbooks on experimental methods 
(see, for example, (2), (3), (4)) and it can there- 
fore be assumed that the basic principles of 
their operations are known. However, it seems 
appropriate to give a summary of the problems 
that arise when the above methods are employed 
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for velocity measurements and to point out their 
major disadvantages : 


1. No local velocity measurements are possible. 
All three methods only allow the determina- 
tion of an integral value over the whole 
velocity field. The signal represents, therefore, 
a mean value over the whole depth of the 
test section. 

2. The flow situations along the side walls of 
test section influence the signal in a way that 
is very difficult to quantitatively assess. 

3. The velocity is measured indirectly. Only 
those flow problems that result in density 
changes, and hence refractive index changes, 
can be investigated since only these will 
create the required bending of the light or 
the changes in wave length, resulting in a 
phase shift. 

4. All three methods are suitable to obtain 
qualitative information. They are less suited 
to yield quantitative results. This is partly 
due to the fact that only integral values are 
obtainable and that these values have to be 
obtained from photographic records. 

5. The application of all three methods is so far 
limited to laboratory use only. This is partly 
due to the difficulty in obtaining quantitative 
information and partly due to the fact that 
all optical arrangements require a conside- 
rable amount of space. 


The most important advantage of any optical 

velocity measuring technique is that, ideally, no 
measuring probe has to be introduced into the 
flow and, therefore, no flow distrubances are 
created. This is the main reason why many 
research workers have tried to develop optical 
velocity-meters that do not have the disadvan- 
tages given above. 
Examples of the optical arrangements of such 
devices are shown in Figs. 2— 8. Brief descri- 
ptions of the different devices and their prin- 
ciples of operation are given below. 


VELOCITY-METER OF THOMPSON (5) 


With the velocity-meter of Thompson, see 
Fig. 2, one is actually measuring the time a 
tracer particle needs to pass the space between 
two light marks that can be conveniently created 
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by means of a Koesters prism. The incident 
light of both beams is blocked off by means of 
a mask in front of a light detecting system. 


MIRROR 
PHOTOMULTIPLIER E LASER 
MASK WITH PRA | 
TWO SLOTS TEST SECTION 
/ LENS L 2 


| ,, R 
TE! ISM 
4 LIGHT KOSTERS PRISM 
/ STOP 
LIGHT COLLECTINS 
LENS 


kig. 2 — Optical velocity meter of D. H. Thompson (5) 


A tracer particle scatters light as it passes through 
the two light marks. The scattered light is col- 
lected by a lens and detected by means of a 
photomultiplier in the form of two voltage pulses 
at the photomultiplier output. By measuring the 
time interval between the pulses, and knowing 
the separation of the two laser beams, the par- 
ticle velocity can be calculated. A similar device 
has been proposed by Hausler and Stumm (6). 
In their set-up the time is measured over which 
a certain scattering target is scattering light. 
We can see from Fig 2 that this type of optical 
velocity meter does not measure integral values. 
However, the spatial resolution of the instru- 
ments is too poor for fluid flow research. Fur- 
thermore, velocity measurements employing these 
techniques turn out to be rather laborious. 


PARTICLE VELOCITY-METER OF GASTER (7) 


This instrument measures the velocity of small 
scattering targets suspended in the flowing fluid. 
Light reflected from these targets illuminated by 
a thin sheet of light is focussed onto a trans- 
mission grating. This light passes through the 
grating onto a photomultiplier which is masked 
so that only a restricted volume of fluid is 
sampled. The photo-electric signal generated by 
the passage of an insulated particle image across 
the grating is periodic whith a frequency which 
is equal to the number of grating lines crossed 
in a given time. Thus, for a given optical 
geometry and grating size, the frequency is pro- 
portional to the velocity of the particle. 
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The instrument shown in Fig. 3 employs a 
moving grating. This not only enables the device 
to sense the flow direction but also to control 
the signal frequency. The difference between 
particle and grating velocity is measured by 
means of two counters to yield instantaneous 
velocities. 
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Fig. 3 — Particle velocity meter of H. Gaster 7) 


CROSSED-BEAM CORRELATION TECNIQUE 
OF FISHER AND KRAUSE (8) 


The optical arrangement shown in Fig. 4 has 
been used to study the influence of turbulence 
on the actual mixing process. The method is 
based on light intensity recording of two light 
beams which traverse the entire flow in two 
mutually perpendicular directions. The radiation 
employed is arranged so that it is partially 
absorbed by a constituent of the flow. Thus, 
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Fig. 4 — Crossed-beam correlation technique of M. J. 
Fisher (8) and F. R, Krause 
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turbulence-induced fluctuations of either the 
thermodynamic properties or the concentration 
of the chosen constituent are reflected as 
fluctuations of intensity of the resultant beam. 
Each detector records the integral of all the 
fluctuations along the entire beam path. Howe- 
ver, as Fisher and Krause showed, the correla- 
tions of the two signals yields local information. 
The above derivations have shown that the 
crossed-beam correlation method is not suited 
to velocity measurements, since only scalar 
quantities are influencing the signal. Those 
measurements turn out to be not too laborious. 


LIGHT INTENSITY TECHNIQUE OF ROSENS- 
WEIG, HOTTEL AND WILLIAMS (9) 


This method has also been employed to study 
turbulent mixing. One can see from Fig. 5 that 
with this technique the light intensity which is 


Fig. 5 — Light intensity technique of R,E Rosensweig, 
H. G. Hottel, G. C. Williams (9) 


recorded is scattered from a fixed volume into 
the direction of the photo-detector. This light 
intensity is ideally proportional to the concentra- 
tion of scattering particles within the control 
volume. Thus, fluctuations of concentrations of 
scatterers are recorded as fluctuations of light 
intensities. 

This light intensity technique is again not 
suited to velocity measurements but represents 
a very useful method to measure concentration 
fluctuations. 


FRINGE ANEMOMETERS 


This group of optical velocity-meters has often 
been referred to as double-beam Doppler-shift 
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anemometers. The author prefers, however, to 
denote these optical velocity-meters, three exam- 
ples are shown in Figs. 6, 7 and 8, as fringe 
anemometers, since the actual rate is measured 
at which tracer particles cross a set of fringes. 
The distinction between fringe anemometers and 
Doppler-shift anemometers is based on conside- 
rations given in (12). 
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Fig. 6 — Fringe anemometer of Mazumder and Wankum (ro) 


The fringe anemometer of Mazumder and 
Wankum (10), see Fig. 6, employs two intense 
light beams that are obtained by a beam splitter 
prism. Both beams are focussed by a single lens. 
At the focal point of this lens, a fringe pattern 
of known geometry is produced that is used for 
velocity measurements. Crossing particles scatter 
light from each fringe within the pattern and 
this light is focussed onto a photo-multiplier. 
The frequency of the voltage variations can be 
used for velocity measurements. 


Fig. 7 — Fringe anemometer of M. J. Rudd (11) 


The fringe anemometer of Rudd (11), see 
Fig. 7, employs a lens and two slots to produce 
a set of fringes at the focal point of the lens, 
This fringe pattern is reproduced on the photo- 
detector by means of a second lens. With this 
optical arrangement the light blocked of by 
small tracer particles is measured, 
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Fig. 8 — Fringe anemometer based on Fresnel biprism 


Another optical arrangement for a fringe 
anemometer is shown in Fig. 8. This arrangement 
is presently under investigation at Imperial 
College. The set of fringes is produced by a 
Fresnet biprism and then focussed down to form 
a line of bright spots. The light scattered by 
tracer particles when passing the light marks is 
recorded to yield the required velocity information. 


BASIC THEORY OF DOPPLER-SHIFT ANE- 
MOME TRY 


The above introduction to several methods of 
optical velocity measurements did not include 
laser Doppler-shift techniques. Such techniques 
have recently been suggested as the basis for 
optical velocity-meters suitable for local velocity 
measurements in transparent gas and liquid flows. 
Test measurements employing laser Doppler- 
-shift anemometers suggest that this new measu- 
ring technique is a very promising method for 
many measurements in fluid flow research. 
However, laser Doppler-shift anemometry is still 
in the developing stage. 

Detailed knowledge of the basic principles of 
the Doppler effect is required to understand 
Doppler-shift anemometry. It seems therefore 
appropriate to derive briefly the equations that 
will be needed for more profound derivations 
at a later stage. 

There are two phenomena that have to be 
distinguished quite clearly: the frequency shift 
due to a moving source and the one due to a 
moving receiver. Both phenomena result in a 
change of frequency. However, they are different 
in nature. This is indicated in Fig. 9. In the 
first case, where the wave source is moving, the 
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Fig. g — Sketch to explain basic principles of the 
Doppler effect 
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stationary observer will record a change in wave 
length, hence in wave frequency. This is due to 
the fact that the waves are squeezed into a space 
given by the distance (c—U, I,).t whereas, for 
a stationary source, the equivalent distance is 
given by c. t. The wave-length for the moving 
source is: 


= ———— (1) 


. . f Cc 
The corresponding frequency is = —: 
yX 


2 
a Tp 


ENA (moving source) (2) 
(Uh) 
c 


The moving observer also records a frequency 
shift if it passes waves emitted from a stationary 
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source. This is due to the observer's movement, 
since it will pass more or less wave trains, 
depending on the direction of its movement, than 
are sent out by the source: 


or, rewritten: 


== vo, (1— aa) 
Cc 


/ 


(moving observer) (3) 


The frequency shifts due to a moving source 
and to a moving observer are the basic principles 
of all the following considerations. 


frequencies received by observer 1 and 2 
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Fig. 10 — Sketch to explain basic principles of laser 
Doppler-shift anemometry 
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The next step in the explanation of the basic 
principles of Doppler-shift anemometry is a sim- 
ple experiment which is outlined in Fig. 10. It 
shows a sender that emits waves of constant 
frequency in a certain direction which is given 
by the unity vector |,. An aeroplane, that con- 
tains a receiver and a sender for the particular 
waves, penetrates the emitted «wave beam». It 
is known from the basic considerations about 
the Doppler effect, that the receiver in the aero- 
plane will record a wave frequency different to 
the one transmitted if the aeroplane has a velo- 
city component in the direction of the unity 
vector |; i. e. a frequency shift is recorded if 
sw 1; Ú, + O. This frequency shift is given by 
equation (3). 


“p =— “e (1 um Ut) (4) 
G 4 


The sender in the aeroplane is now used to 
emit waves in all directions Since the receiver's 
wave has the frequency », the sender can only 
emit waves of this particular frequency. Two 
observers placed at two different positions (A) 
and (B) will record additional frequency diffe- 
rences due to the fact that the sender in the 
flying aeroplane represents a moving wave 
source. The two frequencies can be calculated 
employing equations (2) and (4): 


/ U, À 
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The difference between the two frequencies 
given by (5) and (6) is: 


se (7) 


Equation (7) can now be approximated by the 
following expression if the velocity U, the aero- 
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plane is small in comparison with the wave 
velocity : 


y a é 
pat ou-m-LOs 
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The last equation indicates that velocity measu- 
rements can be made by recording the frequency 
difference of the two waves observed at positions 
(A) and (B). 

Velocity measurements by means of Doppler- 
-Shift anemometers are based on the considera- 
tions given above. The unspecified wave source 
in Fig. 10 is replaced by a light source of suffi- 
ciently high light intensity; thus a laser can be 
conveniently used. This laser emits light waves 
in a certain direction given by the unity vector 
|,, see Fig. 11. 
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Frequency difference between light scattered from a 
particle in two different directions. 


Fig. 11 — The action of scattering particles on the fre- 
quency of laser light waves 
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The aeroplane is replaced by small tracer 
particles that fulfill the action of the receiver 
and the sender in the aeroplane. This is readily 
understood by the following simple considera- 
tions. If small particles penetrate a laser beam 
a certain amount of light is scattered by them. 
Thus, the particles act as «senders». However, 
before that they must have acted as «receivers», 
i. e. their action can be understood as the 
receiving and sending of light waves. 

The above experiment with the aeroplane has 
indicated that velocity measurements by means 
of Doppler-shift techniques require, in addition 
to small tracer particles, a light collecting and 
combining system to eliminate the frequency 
difference between two scattered light beams. 
Additionally, a recording instrument is required 
that is able to record the information contained 
in the resulting optical signal. 

After having decided on the basic optical 
components that are required for velocity measu- 
rements by means of Doppler-shift techniques, 
more detailed considerations are required concer- 
ning the optical signal. It has been shown above 
that the two scattered light waves collected in 
two different directions k, and m, have different 
frequencies, given by equations (5) and (6). The 
difference between the two waves is given by 
equation (8). This frequency difference has to 
be measured to vyield information about the 
velocity component ÚU,n,. The high frequencies 
of the light waves (of the order of 10!º) do not 
allow both frequencies to be measured and then 
subtracted. 

The way in which 4v is usually estimated is 
as follows. Two scattered light waves, y, and y,, 
are superimposed by means of a mirror beam 
splitter arrangement to yield : 


| x 
yY=7, PY «sin(2 mv, t— SE J+ 
" Xo 
+ Ys, Sin (2 7 %t— Et. (9) 


If the light beam of higher light intensity is 
denoted as beam 1, equation (9) can then be 
rewritten as: 


x 
y = [ys,—Ys)] «sin[27(v t— Elk 
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Ya toin [2 Oy t= o) + 
. X5 
«+ sin[27(ht——]) (10) 


This equation can in turn be rewritten (see Durst 
and Whitelaw (12)): 
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The particle acts as light source for scattered waves 
in directions k, and m,. 

x, and x, are changing due to movement of the 
particle: 


x=(0), +0,k;(t—t) 

xo=(x)o +Um,(t—t,) 
(xx) = [xo +] + UK +Am(t—to) 
(x—x9)=[(xDo—(K9)o]+U(k— m)(t—to) 


Fig. 12 — The influence of the change in particle posi- 
tion on the path length of two scattered waves 


Durst and Whitelaw (12) show that equation (11) 
can be written in the following form if the 
considerations indicated in Fig. 12 are taken 
into account. 


y=[y,—ys) -sin[27 (,t—4)] + 
+ 2ys, (sin [27 (1,t— 4,)]) 
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Av A 
[cos [25 (t+ t— 49) (12) 
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this term results from this term results from a 
the Doppler effect change of particle position 


An optical signal as expressed by equation (12) 
will result in an anode current of the following 
type if it is incident on a photo-detector : 


Ys 

2 9 | 
= (Ys—Ys9) + Zn Yo — 

Sg 

— 1] cos[27 (Avt—y,)] + 


+ 2Ys) cos [2x (Avt— Y5)] (13) 


This anode current varies with the frequency 
Av and can therefore be used to measure dv. In 
practice, the frequency of voltage variations 
across a load resistance of a photodetector is 
measured and the resulting signal is processed 
to yield the required information. Four methods of 
signal processing are discussed later in the paper. 


OPTICAL ARRANGEMENTS FOR LASER 
DOPPLER-SHIFT ANEMOMETERS 


Several optical arrangements have been used 
in the past for laser Doppler-shift measurements. 
A short discussion of some typical examples 
might be helpful to demonstrate how several 
research workers have practically solved the 
problems that were indicated during the theore- 
tical treatment of Doppler-shift anemometry in 
the previous section. 

The optical system shown in Fig. 13 has been 
used by Lehmann (13) during his earlier experi- 
ments. With this optical arrangement, scattered 
light is combined with a reference beam that is 
taken from the incident light beam before the 
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Fig 13 — One-dimensional optical arrangement 
of B. Lehamann (13) 
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HERÓIS DA LUBRIFICAÇÃO 


Makerel (Viking) 


Os Vikings foram famosos na pirataria da Iuade Média. 
Mas não se confranjam as consciências puras... À pirataria, 
já nesse tempo, era, além duma lucrativa actividade, um 
título de glória... 

Ole Cavala, chefe de manutenção da frota Viking foi 
o homem responsável pelos êxitos desses intrépidos homens 
do mar. 

As primeiras surtidas dos Vikings tinham sido autênticos 
falhanços e o público só não os assobiou porque não havia 
ainda, como há hoje, público para os Vikings.. 

Os guerreiros remavam nas suas longas embarcações, 
aquelas mesmas que, séculos volvidos, tão popularizadas 
seriam pelo cinema e outros meios audiovisuais de educa- 
cao e diversão. Remavam até aos portos da estranja, pos- 


TÉCNICA-XVII 


did 


à LUBINITILALCAO 


Makerel (Viking) 


sessos da expectativa natural que antecede o fogo-posto, a 
pilhagem e os proventos de um saque em larga escala. 

Infelizmente para eles, as cidades a que arribavam 
tinham sido quase sempre, com prudente antecipação, aban- 
donadas e as presas de guerra (mulheres, na deliciosa maio- 
rio dos casos) há muito que haviam debandado para o inte- 
rior aos gritinhos, mais de excitação que de medo. 

«É por causa da chiadeira que esses malditos toletes (*) 
fazem!», explicou Ole Cavala, o chefe de manutenção da 
frota. «Os habitantes das cidades ouvem-nos a é fiordes 
de distância O nosso problema, é um problema de lubri- 
ficação». 

Ole Cavala sabia que seria preciso consultar um perito, 
para se chegar a uma conclusão definitiva sobre o tipo de 
lubrificação indicado para um caso desses. Por isso, pru- 
dente e avisado como era, resolveu chamar, às boas, um 
fornecedor que, além de vender lubrificantes, também sabia 
resolver problemas com eles relacionados. 

O delegado da MOBIL, contudo, ficou atónito com 
aquele novo tipo de problema. Na sua opinião, nem o 
engenheiro que prestava assistência às vendas poderia resol- 
ver tal bico de obra. E por isso aventou outra hipótese. 

«No nosso laboratório central de investifacão, traba- 
lham inúmeros especialistas de lubrificação», começou ele 
a dizer ao Ole Cavala. «Um deles é perito em toletes rechi- 
nantes e, embora seja um cientista em toda a acepção da 
palavra, é uma pessoa muito simples e acessível. Durante 
três lustros viveu num barco só para estudar toletes. É ele 
o nosso homem». 

Ole Cavala e o delegado delinearam as permissas do 
problema e enviaram-no para a Sede por um falcão-men- 
sageiro. 

Dias depois, o falcão regressou com a resposta: «Expe- 
rimentem Massa Lubrificante de Rena — Mobil N.º 45 Bs». 

O fornecedor, que estava com a mão na massa, pois 
tinha-a em armazém, logo a experimentou com a colabora- 
ção de Ole Cavala. 

E a coisa resultou. Os toletes silenciosos fizeram dos 
Vikings os mais famosos invasores do seu tempo e as mulhe- 
res da estranja passaram a andar num grande desassossego... 

Há muitos, muitos anos já, que Ole Cavala deixou 
de existir, mas a tradição criada por ele ainda vive no espi- 
rito dos chefes de manutenção dos tempos modernos. Eles 
não ignoram que uma parte importante da assistência da 
MOBIL começa nesse grupo de especialistas que trabalham 
nos seus laboratórios. 

Hoje em dia, tais peritos não são muito procurados 
para tirar chiadeiras a toletes, mas a verdade é que sabem 
tudo o que deve saber-se sobre óleos de arrefecimento, 
lubrificantes para motores, rolamentos e outros órgãos de 
máquinas. 

Os vastos e especializados conhecimentos destes 


homens estão à disposição de todos aqueles que compram 
produtos MOBIL. 


Meditem um pouco sobre esta história. É uma 
história que pouca «chiadeira» tem feito, mas 
que se reveste de grande significado! 


(*) Toletes: forquilhas onde se apoiam os remos. 
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measuring point. This results in a decrease of 
light intensity within the measuring control 
volume and, hence, the scattered light intensity 
is also decreased. This is disadvantageous. 
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Fig. 14 — One-dimensional optical arrangement 
of J. W. Foreman et al. (14) 


The optical arrangement sketched in Fig. 14 
has been put forward by Foreman et al. (14). It 
is basically the same as the system of Lehmann. 
However, the reference light beam is taken of 
after the light has passed the scattering region. 
Thus, the total energy of the laser beam is 
employed for the actual scattering process 
resulting in higher signal to noise ratios. 

The instruments shown in Figs. 13 and 14 have 
to be considered as one-dimensional Doppler- 
-shift anemometers, since they are only applicable 
for measurements of a single velocity component. 
Rotation of the whole optical system is required 
if other components have also to be measured. 
This will result in a decrease in the accuracy of 
localisation of the measuring point. 
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Fig. 15 — Two-dimensional optical arrangement ol F, 
Durstand J. H. Whitelaw (1) 


The optical arrangement shown in Fig. 15 has 
been developed by Durst and Whitelaw (1) for 
two-dimensional measurements. This two-dimen- 
sionality is achieved by a simple rotation of the 
mirror-beam-splitter system after the second lens. 
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This rotation does not affect the measuring point 
which is fixed by the light collecting lens L, and 
the two slots. The design of the above systems is 
based on the theoretical considerations given below. 

A three-dimensional Doppler-shift anemometer 
has been developed for NASA by the Raytheon 
Company. The system, as far as the optic is 
concerned, is basically the arrangement used by 
Foreman but three optical set-ups are built in 
three directions. 
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Fig. 16 — Three-dimensional optical arrangement built 
for NASA by Raytheon Company (15) 


The optical system shown in Fig. 17 is again 
a one-dimensional system. The reference beam 
is taken off in a similar way as shown in Fig. 13. 
However, a special combining system is chosen 
to superimpose two light waves. The reference 
beam is passed through the measuring control 
volume and, thus, the light scattered in the 
direction of the reference beam interferes to 
yield the required beat frequency. 


LIGHT SOURCE AND OPICAL GEOMETRY 


The different optical geometries discussed in 
the previous section show that several optical 
systems have been used for laser Doppler-shift 
measurements. This might give rise to questions 
concerning the influence of the geometries on 
the measured frequency. Hence, the question 
might be asked whether there exists an optical 
system. 

The following derivations represent simple 
theoretical considerations concerning the require- 
ments that have to be fulfilled to obtain an 
optimum precision from a laser Doppler-shift 
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Vig. 17 — One-dimensional optical arrangement of Goldstein et al. (16) and Pike et al. (17) 


anemometer. All these derivations are based on 
the final evaluation equation which has been 
derived during previous considerations, see 
equation (8). This equation is valid for all optical 
geometries. However, each optical arrangement 
has its particular »n,. 

It has been common practice to describe velo- 
city components of turbulent flows in terms of 
a mean velocity and a fluctuating component. 


(14) 


If a time-averaging procedure is applied to a 
combination of equations (8) and (14), then the 
following relationship can be derived: 


am 1 
Av == do (U, Dn.) (15) 


This equation implies that, for all optical geo- 
metries, the mean velocity can be calculated by 
measuring a time-mean value of the Doppler 
frequency. Thus, as far as mean velocity measu- 
rements are concerned, all optical systems are 
equivalent. 

One of the great advantages of Doppler-shift 
anemometiy is that it can also be employed to 
measure turbulence characteristics. This can 
readily be seen from the following equation 
which is derived from the evaluation of equa- 
tion (8): 


y 238 1 
o ao a dd U,9n, + — U; n,9» 
c 


/ | É (16) 
| | 


velocity variations chromatic coherence 
optical imperfection 
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The last equation indicates that variations in 
velocity, as they appear in turbulent flows yield 
a variation in frequency. Hence, the spread of 
the measured frequency spectrum can be used 
to obtain turbulence characteristics. However, 
equation (16) shows that this spread is influenced 
not only by velocity variations but also by varia- 
tions of other quantities; e.g. variations of », 
and v. These additional influences, and the 
reasons for their appearance, have to be tho- 
roughly understood to be able to eliminate them 
or to decide whether they can be neglected. Phy- 
sically the two terms express influences due to 
changes or due to imperfections in the light 
collecting system, »,/c. U, ànp;, and due to fre- 


. * - . 1 * 
quency variations in the light source, U,n,%.. 
c 


CHROMATIC COHERENCE 


The relationship between frequency spread and 
variations of different quantities appearing in the 
evaluation equation (8) can be used to derive 
theoretically the basic requeriments for the 
chromatic coherence of light sources employed 
for Doppler-shift measurements, The require- 
ments for chromatic coherence have to be such 
that the spread of the measured frequency spec- 
trum is not influenced by the frequency spread 
of the incident light spectrum. It is sufficient for 
all pratical applications that the last term of 
equation (10) be small in comparision to the 
first one; hence: 


(17) 


The expression on the righthand side of equation 
(17) can be kept small even when conventional 
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ight sources are employed; e.g. characteristic 
values for the green mercury line are: 


(18) 


Thus, Doppler-shift anemometer measurements 
do not require lasers to be used for chromatic 
coherence. The real reason for the employment 
of lasers is their high light intensity which is 
required to obtain measureable scattered light 
signals and their long coherence length. The 
latter makes experimental investigation easier 
since requirements concerning optical path diffe- 
rences are practically eliminated. 


OPTICAL IMPERFEC TION 


The above considerations have shown that, for 
all practical purposes, equation (16) can be writ- 
ten in the following approximated form: 


(Av) = Sn, 3U, + Un) (19) 
c 


rd a 


frequency spread due to frequency spread due to 
valocity variations optical imperfection 


The physical meaning of the two terms of equa- 
tion (19) has to be understood to be able to 
optimise the optical geometry of Doppler-shift 
anemometers: 


1. The first term on the righthand side of 
equation (19) represents the phenomena one 
wants to measure; i.e. the spread of the 
frequency spectrum can be used to obtain 
information about turbulence. 

2. The second term represents an additional 
spread of the measured frequency spectrum, 
This spread should be avoided, since it is not 
due to variation of velocity but to a certain 
imperfection in the optical system. 


Fig. 18 outlines the meaning of the optica] 
imperfection mentioned in 2 above. The light 
collecting system always contains finite aper- 
tures that are equivalent to variations of 
k, and m, in our theoretical treatment of 
Doppler-shift anemometry. This will consequen- 


tly lead to variations of n, and, hence, to a 
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frequency spread that is not due to velocity 
variations. This can lead to considerable influ- 
ences on the signal if the velocity U, is large. 
In this case U, ôn, will be large, and n, 2U, 
may be small. Thus, the influence of the 
optical imperfection tends to plot out the 
spread due to turbulence. In the seconds ketch 
of Fig. 18, the variation of n, has been 


” ” 
+ — "4! mão 


Fig. 18 — Sketch to demonstrate that preferable light 
combining systems exist 


reduced by combining the outer part of the 
upper beam with the inner part of the lower, 
and vice versa. This leads to a smaller value of 
ôn, and, hence, to a reduction of the influence 
of the optics on the signal. 

Fig. 19 shows the pratical solution to the 
problem that is pointed out in Fig. 18. This solu- 
tion is based on the fact that one can control the 
variation of nm. The upper sketch shows an 
optical beam combining system which is based 
on the first beam arrangement show in Fig. 18. 
I will therefore strongly influence the signal, 
ie. U, 9n, is large. The lower sketch of Fig. 19 
shows that, with the chosen mirror-beam-combi- 
ner-arrangement, the outer part of the upper 
beam is combined with the inner part of the 
lower beam, and vice versa. This arrangement 
minimises the variations of n, and, hence, the 
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Optical arrangement leads to increased spread 
of frequency spectrum 
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Optical arrengement with minimised spread 
of frequency spectrum 


Fig. 19 — Good and bad optical light 
combining arrangements 


frequency spread. It should be pointed out that 
the spread is not eliminated but only reduced 
to an acceptable minimum. The minimum of 
frequency spread due to the optics that can be 
accepted is given by the following relationship: 


Un, <<] 1,:9U, (20) 


Í i j Í 


ELECTRONIC DATA PROCESSING 


So far the basic problems indicated have been 
those involved in the design of an optimum 
optical system for laser Doppler-shift anemc- 
metry. However, the optical problems are not 
the only ones that occur. The optical signal has 
to be transferred into an electrical one and this 
signal has to be processed to yield the required 
information. This involves quite a few electronic 
problems that will not be dealt with in this 
paper. However, it seems appropriate to mention 
some methods of data processing that are 
presently being investigated and developed at 
Imperial College. 

All data processing systems have in common 
the fact that a photomultiplier is used to trans- 
form the optical signal into an electrical one, 
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The frequency of the output voltage across the 
load resistance of the photomultiplier contains 
the required information; the Doppler-shift 
frequency. However, this instantaneous infor- 
mation has to be processed to yield the mean 
velocity components and turbulence character- 
istics. 


1. Frequency Analyser 

This method of data processing results basically 
in a frequency probability distribution which 
corresponds tothe velocity probability distribution 
at a particular measuring point. All the required 
information can be obtained from this distribution. 
However, the accuracy is rather limited and the 
evaluation of the data proves to be rather 
laborious. 


2. On-Line Computer 

Small on-line computers can be employed to 
obtain conveniently mean velocity information 
and mean values of the turbulence fluctuations. 
This evaluation method is very accurate. However, 
good signal to noise ratios are required if on-line 
computers are employed, or the signal has to be 
separated before processing it by means of 
computers. 


3. Frequency-to-voltage conversion 

A commercial frequency-to-voltage converter 
can be employed to yield the mean velocity and 
the fluctuating components. However, most of 
them have a very limited linear range. 

New frequency-to-voltage converters have been 
put forward which show a larger linear range. 
These devices can be employed successfully if 
large signal-to-noise ratios can be achieved. 
Otherwise, it must be expected that voltage 
outputs result from the noise and not from the 
signal. 


4. Frequency tracking discriminators 

Low signal-to-noise ratios, as they appear in 
«clean» air measurements, present a big problem 
for data processing of Doppler-shift information. 
The only devices that seem suitable for these 
low signal values, as well as large signals, are 
frequency tracking devices. Such devices eliminate 
most of the disturbing noise by means of filters 
and a kind of loop system is used to pick up 
the varying signal. 
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Fig. 20 — Characteristics of particles and particle dispersoids 


TECNICA N.º 408 443 


THE SIZE OF SUITABLE TRACER PARTICLES 


It has been pointed out during the theoretical 
treatment given above, that laser Doppler-shift 
measurements are based on the scattering action 
of small tracer particles. Therefore, the accuracy 
of measurements of mean velocity and of turbu- 
lence characteristics depends on the accuracy 
with which the velocities of inserted particles 
represent, at any instant, the velocity of the 
flowing Ffluit. Thus, there is a need for tracer 
particles of a carefully predetermined size. 
Whether these particles are already in the fluid 
or have to be added, depends on the nature of 
the fluid chosen for the experiments. 

In practice there always exists a certain particle 
range. The lower limit of this size range is given by 
the fact that small particles are subject to 
Brownian motion and represent therefore the flow 
velocity and a superimposed velocity spectrum due 
to the thermal random motion of the molecules 
and/or atoms. The rms-value of the tracer particle 
movement increases with decreasing particle size. 
The upper limit of the particle size range is given by 
the inertia of the particles which increases with 
increasing particle size. The size of tracer partic- 
les, already in the fluid or artificially introduced, 
has to lie within the above range. 

Durst and Whitelaw (1) derived the limits of 
the appropriate particle size range for turbulence 
measurements in gas flows. They concluded from 
available investigations of particle movement, 
that their particles should be of the following 
size : 

0.1 gm <r<lIium 


In Fig. 20 some characteristics of particles and 
particle dispersions are given to assist in 
choosing the right means of particle production 


FINAL REMARKS 


I seems appropriate to summarise the major 
advantages of laser Doppler-shift tehcniques for 
local velocity measurements in transparent fluids. 


a) Laser Doppler-shift anemometry can be em- 
ployed for both mean velocity estimations and 
measurements of turbulence characteristics. 

b) Instruments based on Doppler-shift techni- 
ques are stable since the only values that are 
in the equation of evaluation are the wave 
length of the incident light beam and the 
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size of the angle between the two directions 
of light collection, 

c) Very localised measurements are possible (the 
dimensions of the measuring control volume 
are of some microns), although improved 
optical systems are required to obtain good 
localisation also for more dimensional mea- 
surements. 

d) The theory necessary for the evaluation of 


experimental data has been derived and has 
proved to be relatively simple. The derivations 


suggest, however, that the processing of 
large quantities of data will require advanced 
data processing systems. 

e) A minimum disturbance of the flow system 
is caused, since the inserted «probes» are 
small particles, e.g. 01 — 1 um, for measure- 
ments in air. 

f) If for any special flow situation the particle 
size has been estimated and if the particle 
material was chosen such that influences of 
temperature, chemical reaction, evaporation, 
etc. can be neglected, the turbulence measure- 
ments by Doppler-shift anemometers are 
independent of temperature and density 
gradients. The influence of velocity gradients 
can be made negligible by collecting light 
from small measuring control volumes. 

g) Measurements in flow configurations of extre- 
me temperature environments are possible. 

h) The inertia of the particles can be kept small 
and the response of the recording equipment 
is practically instantaneous to even the most 
rapid fluctuations. 

The main disadvantage of the laser anemometer 

is that it is still in the development stage. 
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NOMENCLATURE 


velocity of light 

unity vector from particle to observer 
unity vector in the direction of incident 
light beam 

unity vector from particle to observer 
vector difference (k, —m,) 

radius of scattering particles 

time 

instantaneous velocity vector 

mean velocity vector 

fluctuating component of velocity vector 
distances from scattering particle to 
light detector 

amplitudes of light waves 

amplitudes of scattered light waves 
frequency difference (»k, — vm,) 

time mean value of 4» 

variation in 4» (frequency spread) 
variations in velocity (turbulence) 


variation in light frequency of source 
(cromatic incoherence) 

isotropic third order tensor 

wave length of light or of waves 
wave length of Doppler-shifted waves 
wave length of incident light waves 
frequency of light or of waves 
frequency of Doppler-shifted wave 
trequency of light scattered in direc- 
tion k, 

frequency of light scattered in direc- 
tion m, 

frequency of incident light waves 
frequency as «recorded» by moving 
object, e. g. plane, particle 
frequency of light wave 1 
frequency of light wave 2 
factor in equation (13) 
phases of light waves 
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